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We study the quantum-mechanical generation of gravitational waves during inflation on a brane 
embedded in a five-dimensional anti-de Sitter bulk. To make the problem well-posed, we consider 
the setup in which both initial and final phases are given by a de Sitter brane with different values 
of the Hubble expansion rate. Assuming that the quantum state is in a de Sitter invariant vacuum 
in the initial de Sitter phase, we numerically evaluate the amplitude of quantum fluctuations of the 
growing solution of the zero mode in the final de Sitter phase. We find that the vacuum fluctuations 
of the initial Kaluza-Klein gravitons as well as of the zero mode gravitons contribute to the final 
amplitude of the zero mode on small scales, and the power spectrum is quite well approximated by 
what we call the rescaled spectrum, which is obtained by rescaling the standard four-dimensional 
calculation following a simple mapping rule. Our results confirm the speculation raised in Ref. [2^| 
before. 

PACS numbers: 04.50.+h, ll.10.Kk, 98.80.Cq 



I. INTRODUCTION 

In recent years there has been a growing attention to 
braneworld scenarios, in which our observable universe 
is realized as a brane embedded in a higher dimensional 
bulk spacetime p| . Among various models a simple setup 
called the Randall-Sundrum type II model is of partic- 
ular interest because, despite gravity propagating into an 
infinite extra dimension, four-dimensional general rela- 
tivity can be recovered on the brane due to the warped ge- 
ometry of anti-de Sitter spacetime dUQ. Cosmological 
consequences of the Randall-Sundrum type braneworld 
have been widely discussed la, m B. E3 , especially 
regarding cosmological perturbations E2 . 

The generation and evolution of perturbations are 
among the most important issues in cosmology because 
of their direct link to cosmological observations such as 
the stochastic gravitational wave background and the 
temperature anisotropy of the cosmic microwave back- 
ground (CMB), by which we can probe the early uni- 
verse. While the cosmological perturbation theory in the 
conventional four-dimensional universe is rather estab- 
lished [l2L IT^ | , calculating cosmological perturbations in 
the braneworld still remains to be a difficult problem. Al- 
though some attempts have been made concerning scalar 
perturbations El E O El 13, further progress 
is awaited to give a clear prediction about the CMB 
anisotropy in the Randall-Sundrum braneworld. Almost 
the same is true for g ravitational wave (tensor) pertur- 
bations @, 123, H3, H^H^, H3, H3, IIH Ha 113, Hi| • 

However, since their generation and evolution depend ba- 
sically only on the background geometry, they are slightly 
easier to handle. 

During inflation super-horizon gravitational wave per- 
turbations are generated from vacuum fluctuations of 
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gravitons. The pure de Sitter braneworld is the spe- 
cial case that allows definite analytical computation of 
quantum fluctuations. Thanks to the symmetry of the 
de Sitter group, the perturbation equation becomes sep- 
arable and can be solved exactly [2(j. The time varia- 
tion of the Hubble parameter generally causes mixing of 
a massless zero mode and massive Kaluza-Klein modes, 
and this effect was investigated based on the "junction" 
models, which assume an instantaneous transition from 
a de Sitter to a Minkowski brane [2l| or to another de 
Sitter brane p^ |. These works discussed the quantum- 
mechanical generation of gravitational waves within such 
limited and simplified models. 

As for the classical evolution of gravitational waves 
during the radiation dominated epoch, several studies 
have been done [2! 0, HI El , assuming that a given 
single initial mode for each comoving wave number domi- 
nates. The focus of these works is mainly on the evolution 
of modes which re-enter the horizon in the high-energy 
regime. While the late time evolution is worked out an- 
alytically in Refs. 0,03 by resorting to low-energy ap- 
proximation methods. 

In this paper we consider the generation of primor- 
dial gravitational waves during inflation in more general 
models of the Randall-Sundrum type. For definiteness, 
we adopt a simple setup in which both initial and fi- 
nal phases are described by de Sitter braneworlds. Two 
de Sitter phases with different values of the Hubble ex- 
pansion rate are smoothly interpolated. This work is an 
extension of the work on the "junction" models by the 
present authors in collaboration with H. Kudoh [23. In 
the previous analysis the transition of the Hubble rate 
was abrupt and the gap was assumed to be infinitesi- 
mal. Here we extend the previous results to more gen- 
eral models with smooth transition by using numerical 
calculations with a refined formulation. 

This paper is organized as follows. In the next section 
we briefly summarize past studies |20t l2ll |22j on the gen- 
eration of gravitational waves via quantum fluctuations 
during inflation in the Randall-Sundrum braneworld, em- 



2 



phasizing the mapping formula introduced in Ref. |22j . 
In Sec. 11111 we describe our numerical scheme to investi- 
gate the generation of gravitational waves, and then in 
Sec. HVI we present results of our calculations. Section Ivl 
is devoted to discussion. 



II. GRAVITATIONAL WAVES IN 
INFLATIONARY BRANEWORLD 

A. Pure de Sitter brane 

The background spacetime that we consider is com- 
posed of a five-dimensional AdS bulk, whose metric is 
given in the Poincare coordinates by 



ds 2 



-dt 2 



5ijdx z dx-' + dz 2 



(1) 



and a Friedmann brane at z = z(t). Here £ is the bulk 
curvature length, which is constrained by table-top ex- 
periments as i < 0.1 mm |32| . 

First let us consider the generation of gravitational 
waves from pure de Sitter inflation on the brane |20j| . 
The coordinate system appropriate for the present situ- 
ation is 



ds 2 = 



sinh £ 



- 2 {-dr? 
r\ l v 



Sijdx^x 3 ) + d£ 2 



(2) 



which is obtained from Eq. by a coordinate transfor- 
mation 



t = 7/ cosh £ + to, 
z = —r) sinh£, 



(3) 
(4) 



where io is an arbitrary constant and rj is the conformal 
time, which is negative. The de Sitter brane is located 
at 

£ = £b = constant, (5) 

and the Hubble parameter on the brane is given by 

F = ^ 1 sinh£ b . (6) 

The gravitational wave perturbations are described by 
the metric 



ds" = 



1 



[-dry 2 + (fy + h^Ul.r'd.r' + df \ ,(7) 



sinh £ [rj 



and we decompose the perturbations into the spatial 
Fourier modes as 



V2 



(27tM 5 ) 3 / 2 



d 3 M k (r/,£K k - x e, 



(8) 



where ey is the transverse-traceless polarization tensor 
and M5 is the fundamental mass scale which is related to 
the four-dimensional Planck mass Mpi by £(M^) 3 = Mp V 



The prefactor \[2j (M5) 3 / 2 is chosen so that the kinetic 
term for 0k in the effective action is canonically normal- 
ized. Hereafter we will suppress the subscript k. The 
linearized Einstein equations give the following Klein- 
Gordon-type equation for 0: 



d 2 



2d ,0 sinnV d 
-r^ + kr- 



1 



d 



dr/ 1 rj drj 



V 2 d£ sinh 3 e d£ 



= 0, (9) 



which is obviously separable. Assuming the Z2- 
symmetry across the brane, the boundary condition on 
the brane is given by 



n M 0| b ran e OC 8, 



0. 



(10) 



where n M is an unit normal to the brane. Equation @ 
admits one discrete zero mode, 



H 



-ikrj 



(ii) 



where Af is a normalization constant, which will be fixed 
later, as well as massive Kaluza-Klein (KK) modes, 

where ip u (rf) and Xu(€) obey 

' d 2 2d l2 v 2 + 9/4 

L h z J ' 

dr/ 2 



sinh 3 £ 



2d l2 v 2 
r\ dr\ 
d 1 d 



Mv) = o, (12) 
x„(0 = o, (13) 



'd£sinh 3 £c>£ 

and the corresponding Kaluza-Klein mass is given by 
1 9 \ __9 9 



m 2 = 



H 2 > 



:H 2 . 



(14) 



Due to the mass gap Am = 3H/2 @ between the zero 
mode and the lowest KK mode, all the massive KK modes 
decay on the brane at super-horizon scales during infla- 
tion. The explicit forms of the KK mode functions are 
shown in Appendix IaI 

In inflationary cosmology, fluctuations in the graviton 
field (and other fields such as the inflaton) are consid- 
ered to be generated quantum-mechanically. Following 
the standard canonical quantization scheme we quantize 
the graviton field. Treating as an operator, it may be 
expanded as 

/>oo 

4> = ao^o + + dv (a v (f> v + a*0*) . (15) 



The annihilation and creation operators a n and a) n with 
77 = or v satisfy the usual commutation relations and 
the modes 0o and <p v are normalized so that they satisfy 



Wo ■ <Po) = 
{<t>v - 

(00 • <S>v) = 
(4>n ■ 0*') 



-($> ' $>) = 1) 
-(#•#,) = *(!/■ 

(0S • <t>V) = 0, 

0, for 77, 77' = 



(16) 



0,i/, 
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where ( • ) is the Wronskian defined by |2lJ, |2 



(X ■ Y) := -2i [ 



2 ■ , 3 t ( xd n Y * ~ Y*d„X) 



(17) 



From the conditions l|16|) we can determine the normal- 
ization of the zero mode as 



Af = C{£H), 



(18) 



with 



C(x) 



\j\ + x 2 + x 2 In 



1 + s/l + x 2 



1-1/2 



, (19) 



which is the same as that introduced in Refs. |20ll21ll22| . 
and behaves like C 2 (x) ~ 1 for x -C 1 and C 2 (x) « 
3x/2 for a; 3> 1. This function plays an important roll 
throughout the present paper. 

The expectation value of the squared amplitude of the 
vacuum fluctuation in the zero mode is given by 



b \ 2 = £-i C 2 (£H) — (l + k 2 r, 2 ) 



H 2 

; 2fc 3 ' 



(20) 



where the expression in the second line is obtained 
by evaluating the perturbation in the super-horizon 
regime/at a late time, and hence this is the amplitude 
of the growing mode. In terms of the power spectrum 
defined by 



V 



4nk 3 



(2tt)3 (M 5 )3' 



we have 



v= 2C 2 (£H) (HV 



M 2 
JW P1 



(21) 



(22) 



where we have used ^(Ms) 3 = Mpy. This is the stan- 
dard flat spectrum up to the overall factor C 2 (£H). In 
conventional four-dimensional cosmology the power spec- 
trum of the primordi al g ravitational waves from de Sitter 
inflation is given by [33 



V 



4D 



(23) 



Thus just by rescaling the amplitude in four-dimensional 
cosmology as 



H i ► HC{£H), 



(24) 



the braneworld result (1221 is exactly obtained 



B. The "junction" model 

Pure de Sitter inflation on the brane described in the 
previous subsection is a special case where the amplitude 
of the growing zero mode can be obtained completely 
analytically. As a next step to understand gravitational 
waves from more general inflation with H ^ const, a dis- 
continuous change in the Hubble parameter was consid- 
ered in Ref. In such a "junction" model the Hubble 
parameter is given by 



H(t)) = 



and 



Hi, 

H f =Hi- 6H, 



SH i 

H <<h 



V < Vo, 

V > Vo, 



(25) 



(26) 



is assumed. 

The evolution of gravitational wave perturbations can 
be analyzed by solving the equation of motion back- 
ward using the advanced Green's function [2l]]. A set 
of mode functions can be constructed for the initial de 
Sitter stage, and another for the final de Sitter stage, ei- 
ther of which forms a complete orthonormal basis for the 
graviton wave function. The final (growing) zero mode 
may be written as a linear combination of the initial zero 
and KK modes. The Bogoliubov coefficients give the 
creation rate of final zero mode gravitons from the initial 
vacuum fluctuations in the Kaluza-Klein modes as well 
as in the zero mode. Thus the power spectrum may be 
written as a sum of two separate contributions, 



V = Vo + V 



KK, 



(27) 



where Vo and Vkk are the parts coming from the initial 
zero and Kaluza-Klein modes, respectively. 

An important result brought by analyzing this junc- 
tion model is that V is well approximated by the rescaled 
spectrum V ICS obtained by using a simple map speculated 
by the exact result of pure de Sitter inflation [Eq. (|24() ]. 
More precisely, let V^ik) be the gravitational wave spec- 
trum evaluated in the standard four-dimensional infla- 
tionary universe with the same time dependence of the 
Hubble parameter l(25fr. and then the rescaled spectrum 
is defined by 



4D 



M 2 
JU P1 



V re 



2C 2 (£h k ) 



(28) 



Comparing the rescaled spectrum with V obtained from 
a five-dimensional calculation, one finds that the differ- 
ence between these two is suppressed to be second order 
like El 



V - V Te 



V 



< 



SH 



< I. 



(29) 



It might be worth noting here that the KK contribution 
Vkk is necessary for realizing this interesting agreement 
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between the braneworld result and the (basically) four- 
dimensional result especially at high energies £H ^> 1. 
The above agreement Ij29(l raised a speculation that the 
mapping formula hk > C(£hk) may work with good ac- 
curacy in more general inflation models with a smoothly 
changing expansion rate. The central purpose of the 
present paper is to check whether this speculation is cor- 
rect or not. 



III. FORMULATION 
A. Basic equations 

Now let us explain the formulation that we use to study 
the generation of gravitational waves without assuming a 
pure de Sitter brane. Our formulation is based on double 
null coordinates, which are presumably the most conve- 
nient for numerical calculations. 

The metric Q can be rewritten by using double null 
coordinates 



u = t — z, 
v = t + z, 



in the form of 



ds 2 = 



4/ 2 



(v — u) 2 



-dudv + 5ijdx l dx J ^ 



(30) 
(31) 



(32) 



The trajectory of the brane can be specified arbitrarily 
by 



v = q(u). 

By a further coordinate transformation 

U = u, 
q(V) = v, 

we obtain 

ds 2 = 



(33) 



(34) 
(35) 



U 2 



W) - u] 



[-q'(V)dUdV + S ij dx i dx j ], (36) 



where a prime denotes differentiation with respect to the 
argument. Now in the new coordinates the position of 
the brane is simply given by 



U = V. 



(37) 



We will use this coordinate system for actual numerical 
calculations. 

The induced metric on the brane is 



U 2 

da 2 - 

" ~ [q(V) - V} 2 



[-q'(V)dV 2 + S.jdx'dx 3 ] , (38) 



from which we can read off the conformal time r\ and the 
scale factor a, respectively, as 



dr, = vW) dV, 
21 



a = 



q{V)-V 



(39) 
(40) 



and hence the Hubble parameter on the brane is written 
as 



£H = 

or equivalcntly 
q'(V) 



1 



2vW) 



(41) 



(42) 



Given the Hubble parameter as a function of r], one can 
integrate Eqs. I|39|l and l|42|l to obtain q as a function of 
V. 

When the brane undergoes pure de Sitter inflation (and 
thus q'(V) — constant), the following relation between 
(U, V) and (?7,£) will be useful: 



^ b+ 2 ln U3T/ 



n 



--e-^[(t a -U)(t Q -V)}^ 2 , 



and 



q(V) = e- 2 ^(V -t a )+t, 



(43) 
(44) 

(45) 



The Klein- Gordon- type equation for a gravitational 
wave perturbation <f> in the (U, V) coordinates reduces 
to 



G 



q(V) - U 



(dy — q'(V) djj) + q'(V) k 2 



supplemented by the boundary condition 

[du-dv]4>L^r=0. 



\u=v 



= 0, 
(46) 

(47) 



The expression for the Wronskian evaluated on a con- 
stant V hypersurface is given by 



(X-Y) = 2i dU 



2( 



1 3 



q(V) - U 



(XduY*~Y*duX), 



(48) 



which is independent of the choice of the hypersurface. 

When inflation on the brane deviates from pure de Sit- 
ter one, the decomposition into the zero mode and KK 
modes becomes rather ambiguous. For this reason we re- 
quire the initial and final phases of inflation to be pure 
de Sitter, though arbitrary cosmic expansion is allowed 
in the intermediate stage. In both de Sitter phases, q(V) 
can be fit by Eq. I|45|l . and & and to are determined, 
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respectively. Hence we have two sets of de Sitter coor- 
dinates (?7, £) and (??,£). We distinguish the coordinates 
in the hnal phase by associating them with tilde. In the 
final de Sitter phase the mode will be well outside the 
horizon, and hence we expand the graviton field in terms 
of the growing and decaying zero mode solutions 4>g and 



as 



<f> = Ag<j)g + A d (f> d + \ d,V (a v 4> v + 4 , (49) 



where 



b g := Im 4> 
l d := Re 4> 



(50) 
(51) 




and the mode functions with tilde are defined in the same 
way as 4>$ and <f> v with the substitution of (£ , fj) for (£, rj). 
It can be easily seen that the growing and the decaying 
modes are normalized as 

(^•^) / = ^ (^•0 9 ) / = (^-^) / =O, (52) 

where subscript / means that the expression is to be 
evaluated in the final de Sitter phase. Notice that de 
Sitter mode functions thus defined as functions of (U, V) 
through (£,??) (or (£,?7)) do not satisfy the equation of 
motion outside the initial (or final) de Sitter phase. Back 
in the initial de Sitter phase the graviton field can be 
expanded as 



<j) = a <po + ap^o + / dp (a v <j) v + aj,0*) 



(53) 



We assume that initially the gravitons are in the de Sitter 
invariant vacuum state annihilated by a and a„, 



a |Q) = a„|0) = 0. 



(54) 



We would like to evaluate the expectation value of the 
squared amplitude of the vacuum fluctuation in the grow- 
ing mode at a late time, 



(0|^|0>, 



or equivalently, the power spectrum, 



FIG. 1: Numerical (backward) evolution scheme. 



Suppose that a solution $({/, V) is chosen so as to sat- 
isfy $ = 4>d in the final de Sitter phase. From the Wron- 
skian condition \\W2\\ we see that 



o • $)i a + / dv{4> v ■ $)i a„ + h.c. 



-A n 



where subscript i means that the expression is to be eval- 
uated in the initial de Sitter phase. In the above we used 
the fact that the Wronskian is constant in time. Thus we 
obtain 



(0|A 2 I0) 



1(00-$) 



(56) 



from which we can calculate the power spectrum of the 
primordial gravitational waves. It is obvious that the 
spectrum is written in the form of Eq. I|27|) with 



8C\m f ) (Hj_ 
2tt 



Ml, 



(57) 



V{k) 



Airk 3 



(2tt) 3 (M 5 ) 3 



(0|i 2 J0) 



2C 2 (eH f ) (H f \ /nl u 



Ml, 



2^) (01^10), 



(55) 



where Hf is the Hubble parameter in the final de Sitter 
phase. In order to obtain the final amplitude, it is not 
necessary to solve the evolution of all the (infinite number 
of) degrees of freedom with their initial conditions set by 
Eq. (|53|l . In fact, we have only to solve the backward 
evolution of the final decaying mode, as explained below. 



B. Numerical scheme 

Here we describe the algorithm that we employ to solve 
the backward evolution of the gravitational perturbations 
numerically. We set the boundary conditions for the 
mode $ so as to be identical to a decaying zero mode 
4>d at a time V = Vf in the final de Sitter phase. We 
decompose $ as 



$(U,V) = <f> d (U,V)+6${U,V), 



(59) 



6 



and solve the equation for 6$(U, V) instead of <&(£/, V). 
The equation of motion for <5<£> is obtained as 



Adudy + 



6 



q(V) - U 



(8v -q'{V)du) + q'{V)k< 



5$ 



Adud 



6 



uuv 



q(V) - U 



(dv-q'W^+q'iV) k 2 



(60) 



Since both $ and <pd satisfy the boundary condition of 
the form of Eq. I|47|) . the boundary condition for 5$ is 
also written as 



[da-dv]5Q\ 



u=v 



0. 



(61) 



We immediately find that we do not have to solve the 
backward evolution in the final de Sitter phase, since (5<I> 
identically vanishes there. 

In some cases <E> is not disturbed so much from its fi- 
nal configuration <j)d(U,Vf), especially when we discuss 
a small deviation from the pure de Sitter case. It is ad- 
vantageous then to use a small quantity (5$ as a variable 
in numerical calculations. Of course, there is no problem 
even when does not stay small. 

Our scheme to obtain the values of 5<fr at the initial 
surface V — Vi is as follows. We give the boundary con- 
ditions at V = Vf as /2,i, • ■ • , /jv,i = 0, where 

f n>m := 6<£(U = V f - e(n -l),V = V f - e(m - 1)). 

A sketch of the numerical grids is shown in Fig. 2] Here 
N is taken to be sufficiently large, and we use the same 
grid spacing e both in the U and V directions. At a 
virtual site (U, V) = (Vf, Vf — e) we set 



jfl.2 — II 



(62) 



so that the boundary condition (fol"J) is satisfied at (Vf — 
e/2, Vf — e/2). From fi : i, / 2 ,i, and /i i2 we can determine 
the value of / 2 ,2 by using the equation of motion (|6T7|l at 
(Vf - e/2, Vf - e/2), and then from / 2 ,i, / 3 ,i, and / 2 , 2 
we can determine /3 )2 i and so on. We repeat the same 
procedure in the subsequent time steps until we obtain 
Jm,m, Jm+i,m, ■ ■ ■ , Jn,m at Vi = Vf — e(M — 1). 



C. A toy model 

As stated above, to make the problem well posed, we 
consider models which have initial and final de Sitter 
phases. To perform numerical calculations, as a concrete 
example, we adopt a simple toy model in which the Hub- 
ble parameter is given by an analytic form 



H(rf) = H - Atanh 



where 



H 



Hi + H f 



V~Vo 



Hj — H 



2 ' 



(63) 



(64) 




-100 -80 -60 -40 -20 




FIG. 2: The Hubble parameter and the slow-roll parameter e 
in our model. This is a plot for the model of Fig. 



with Hi and Hf the initial and final values of the Hubble 
parameter, s is a parameter that controls the smooth- 
ness of the transition, and t)q indicates a transition time. 
Taking the initial time r\i and the final time r/f so that 
Vo — Vi 3> s and rjj — rj ps —770 » s, we have H w H t 
for 77 — ► rji and H « Hf for 77 > 77 y . The scale factor for 
this inflation model is given by 



a(rf) = < sA In 



2 cosh 



-Hr] + A770 



(65) 



where the integration constant is determined by impos- 
ing that a sa (—Hf-q)^ 1 at 77 w r)f w 0. The slow-roll 
parameter, e := —d v H/(aH 2 ), is obtained as 



A{sA In [2 cosh((77 - t? )/s)] - Hr] + At] } 
s[H cosh((?7 — t]q)/s) — Asinh((r] — ?7o)/s)] 2 



(66) 



which takes maximum at 77 = r/o- The maximum value is 



sH 2 a(r] ) 



(67) 



The behavior of the Hubble parameter and the slow-roll 
parameter is shown in Fig. [3 

Here we should mention the limitation of this simple 
toy model. For fixed values of Hi, Hf, and e ma x one can 
prolong the period of the transition ~ s measured in con- 
formal time by taking large s. However, using the defini- 
tion of e max (|67l) , the transition time scale in proper time 
is found to be given by ~ a(rjo)s = H 2 e max / A. Therefore 
we cannot change the transition time scale independently 
of the other parameters, Hi, Hf, and e max . 
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FIG. 3: (color online). Top panel: power spectra of gravi- 
tational waves normalized by (2C 2 (ffli) / , M£ 1 ){H i /2n) 2 . Red 
diamonds (upper ones) indicate the result including the con- 
tributions from the initial Kaluza-Klein modes, while green 
diamonds (lower ones) represent the contribution from the ini- 
tial zero mode. Although the rescaled four-dimensional spec- 
trum is shown by blue diamonds, they are almost hidden by 
the red ones since the result of the five-dimensional calculation 
is well approximated by the rescaled four-dimensional spec- 
trum. Bottom panel: difference between the five-dimensional 
and rescaled four-dimensional spectra. 



TABLE I: Parameters used for the numerical calculations 
presented in the figures. 



FIG. 4: Same as Fig. [2] but the parameters are different. 
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FIG. 5: Same as Fig. [2] but the parameters are different. 



IV. NUMERICAL RESULTS 

Using various parameters shown in Table. [I] we per- 
formed numerical calculations, the results of which are 
presented in Figs. 1514101 Numbers of grids are N X M = 
50000 x 1000, and the grid separation e is chosen to be 
(0.11 — 2.5) x £ depending on the energy scale of infla- 
tion. Note here that the step width in conformal time 



rj is given by A?? = y/q'(V) e [Eq. ffity]. and q'(V) be- 
comes smaller for a larger value of £H [Eq. i|42[l ]. Our 
choice of e makes Arj about the same size, Ar/ rs 0.1, in 
all the calculations. Integration over v is performed up 
to v = 20 - 40 with an equal grid spacing of 0.15. The 
typical behavior of the integrand \ {<j) v ■ ^)i\ 2 is shown in 
Fig. EH 
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FIG. 6: Same as Fig. but the parameters are different. FIG. 8: Same as Fig. [3] but the parameters are different. 
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In order to compare the five-dimensional power spec- 
trum with a four-dimensional counterpart, we solve the 
conventional evolution equation for gravitational waves 
in four-dimensions, 

(8g + 2aHd v + k 2 ) $ = 0, (68) 

where H(rj) is given by the same function as used for the 



corresponding five-dimensional computation, and then 
we calculate the rescaled power spectrum obtained from 
the four-dimensional 'bare' spectrum by using the map- 
ping formula l|28l) . 

We introduce a parameter that represents a difference 
between the five-dimensional power spectrum and the 
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rescaled four-dimensional spectrum: 
S( k ) := V ™~ V . 



(69) 



In all the cases of our calculations we clearly see that 
\5\ -C 1 and thus we conclude that the primordial spec- 
trum of the gravitational waves in the braneworld is quite 
well approximated by the rescaled four- dimensional spec- 
trum. We also find that the difference becomes larger for 
larger values of the slow-roll parameter e, but \5\ is no 
greater than £>(10" 2 ) even for models with e max ~ 0.75. 
Notice that the universe is not inflating any more if e is 
greater than unity. At high energies £H > 1 it can be 
seen that dependence of S on £H is weak, while at low 
energies like £H ~ 0.1, 5 is further suppressed compared 
to the high energy cases. 

Contributions from the initial KK modes can be sig- 
nificant for small wavelength modes. For example, for 
the modes with k/a(rjo)H ~ 10 we see that Vkk/V ~ 



0(1O -2 ) at low energies, but the KK contribution can 
become as large as Vkk/V ~ 0.6 in our most 'violent' 
model with £H ~ 10 and e max ~ 0.75. 



V. DISCUSSION 

In this paper we have investigated the generation of 
gravitational waves during inflation on a brane and com- 
puted the primordial spectrum. Extending the previous 
work [23, we have considered a model composed of the 
initial and final de Sitter stages, and the transition region 
connecting them smoothly. We have numerically solved 
the backward evolution of the final decaying mode to ob- 
tain the amplitude of the growing zero mode at a late 
time by making use of the Wronskian method, 

We found that the power spectrum V is well approx- 
imated by the rescaled spectrum V rcs basically calcu- 
lated in standard four-dimensional inflationary cosmol- 
ogy. Here the rescaling formula for the amplitude is 
given by a simple map h% > hkC(ihk), where C is the 
normalization factor of the zero mode. Although the dif- 
ference of the two spectra, 5 :— (V ICS — V)/V, depends 
on the energy scale of inflation and the slow-roll parame- 
ter, our numerical analysis clearly shows that in any case 
the mapping formula works with quite good accuracy, 
yielding \6\ < O(10~ 2 ). This implies that the mapping 
relation between the two spectra holds quite generally in 
the Randall- Sundrum braneworld, which we believe is a 
useful formula. 

We have taken into account the vacuum fluctuations in 
initial Kaluza-Klein modes as well as the zero mode, and 
found that both of them contribute to the final amplitude 
of the zero mode. The amount of the initial KK contribu- 
tion can be large on small scales, and it increases as the 
energy scale of inflation £H becomes higher. This gives 
rise to a quite interesting picture. When the expansion 
rate changes during inflation, zero mode gravitons escape 
into the bulk as KK gravitons, but at the same time 
bulk gravitons come onto the brane to compensate for 
the loss, and these two effects almost cancel each other. 
This seems to happen, irrespective of the energy scale, 
in a wide class of the inflation models even with a not- 
so-small slow-roll parameter. It is suggested that this is 
not the case in the radiation dominated decelerating uni- 
verse H^HHEEll, where the decay into KK gravitons 
reduces the amplitude of the gravitational waves on the 
brane. The junction model of Ref. [2l| joining de Sitter 
and Minkowski branes also show the suppression of the 
gravitational wave amplitude at high energies. 

What is the reason for the remarkable agreement of the 
braneworld spectrum and the rescaled four-dimensional 
spectrum? Extremely long wavelength modes, which 
leave the horizon during the initial de Sitter stage much 
before the Hubble parameter changes, have a squared 
amplitude of (2C 2 (WJ/M^) (Hi / '2n) 2 , and the ampli- 
tude of the perturbations stays constant during the sub- 
sequent stages [3^ | . The same is true for four-dimensional 
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inflationary cosmology, and so the mapping formula is ap- 
plicable to these long wavelength modes. On the other 
hand, a mode whose wavelength is much shorter than the 
Hubble horizon scale at the transition time will feel the 
transition as adiabatic and hence the particle produc- 
tion is exponentially suppressed, (0|AJA 9 |0) » 1, lead- 
ing toP« {2C 2 (£H f )/M<- l )(H f /2ir) 2 . The same argu- 
ment can be applied to the conventional cosmology. Thus 
it is not surprising that the mapping formula works for 
such short wavelength modes. However, at present there 
seems no simple reason for the amplitude of the modes 
with k ~ 0(a(rjo)H) to coincide with the rescaled one. 
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APPENDIX A: KALUZA-KLEIN MODE 
FUNCTIONS 

In the de Sitter braneworld the Kaluza-Klein mode 
functions in the time direction satisfy Eq. i|12fl and are 
given by 

Mv) = ^l-^e-^i-^H^i-kr,), (Al) 



where H^} is the Hunkcl function. They are normalized 
so that 



i£ 3 



(A2) 



The mode functions in the extra direction satisfies 
Eq. i|13|) and are given by 



Xu(0 = Ci(sinh0 2 PrJ^coshO 



C 2 0_i /2+w (coshO 



(A3) 



where P_ 1 / 2+ j I/ and Q_1/ 2 +w are * ne associated Legen- 
dre functions. The constants C\ and C2 are determined 
by the normalization condition 



> (b (sinh£) 3 
and the boundary condition (|47(l . respectively, as 



(A4) 



r(5/2 + if 

T(-iu) 



irC : 



r(«i/-3/2) 



Co = 



T(5/2-iu) " * T(l + iu) 

P -l/2 +l ,( COSh ^) 



-1/2 



-1/2+1 



, (cosh £ b ) 



(A5) 
(A6) 
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